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Observing behaviourally unconstrained animals can lead to 
simple characterization of complex behaviour. We apply 
this principle to infer the neural dynamics of self-initiated 
movement in pulse-type electric fish, Gymnotus sp. 

Recent long-term monitoring of fish (~200 hours from 22 
recording sessions, over 4 animals) in freely swimming con- 
ditions, devoid of external stimuli, reveals non-trivial struc- 
tures in their pattern of electric organ discharge (EOD). 
Simultaneous recording of EODs and fish movement show 
that the EOD rate (EODR) and the activity level of the fish 
are bi-modally distributed, as well as highly correlated. 
These features thus effectively define behavioural attractor 
states corresponding to high and low levels of neural activ- 
ity (up- and down-states, respectively). Trajectories in the 
EODR-activity plane consist of diffusional motion around 
the attractor states, interrupted by sharp transitions 
between states. The duration of each state is uncorrelated 
with that of the next up- or down-state, and is log-normally 
distributed with no characteristic time-scale. 

Based on this data, our goal is to develop a modelling 
framework to better understand the neural pathway 
responsible for self-initiated movement. However, because 
the physiological parameters defining this pathway are 
experimentally unconstrained, it would be premature at 
this point to develop a detailed biophysical model of this 
system. There is thus a preliminary need to, instead, char- 
acterize the key features of the data from a phenomenolo- 
gical perspective. To address this research gap, we attempt 
to fit a stochastic process, with the simplest combination 
of dynamical components, that most closely reproduces 
the statistics of the data. 

As a first step, we hypothesize that the first principal 
component of the data (EODR and activity level) follows 



an overdamped Brownian motion in a double-well poten- 
tial with additive noise. Based on the approach of [1], we 
fit a 4* order polynomial for the potential function by 
associating the stationary solution of the Fokker-Planck 
equation with the measured histogram. For most record- 
ing sessions, the fits appropriately reproduce the histo- 
grams, but consistently underestimates the width of the 
up-state potential well. 

Once the potential function is determined, we generate 
an estimate for the noise intensity by calculating the mean 
escape time from small regions at the bottom of either 
well, and comparing it with its theoretical expression. We 
find, however, that this estimate is dependent on the well 
that was used to generate it, with all recording sessions 
showing a larger noise estimate for the up-state than for 
the down-state. Moreover, in most cases, the ratio of up- 
to down-state duration is underestimated by the fitted 
process, indicative of either too short up-states, or too 
long down-states compared to the data. 

These findings are consistent with the above discre- 
pancy for the width of the up-state potential well, as 
well as with visual inspection of the data, which shows 
greater variability when fish are in up-states. Taken 
together, these observations strongly suggest that state- 
dependent noise is involved in the process generating 
the data, either in the form of multiplicative noise, or 
Poisson shot noise. 
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